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ABSTRACT 


The nuifierical solution of the coupled differential equations of motion of the blades of an horizontal 
axis wind turbine Is a more direct approach than the technique of finite elements, permitting the 
optimization of the design at relatively low cost. The procedure c onsi sts in transforming the equation 
of motion into a set of first order equations and solving them with fourth order Runge-Kutta integrators. 
This technique is applied to a twisted, tapered blade of variable cross section and stiffness including 
discontinuities. The first six natural frequencies and mode shapes are obtained. 

This technique is extended to obtain the polar moment of inertia of the blades as a function of 
frequency and rotational speed. 

A good niatch with the experimental results is achieved. 


INTRODUCTION 

The Accurate determination of natural frequencies 
is of fundamental Importance In the design of 
wind turbine blades. Similarly, the polar moment 
of inertia of the blades is required for the 
study of the torsional dynamics of the drive 
train. 

Rotor dynamics Is often studied with the use of 
large and specialized finite elements computer 
codes. However, the availability and cost of 
operation of these programs limit their use, and 

direct approach could be beneficial* The 

direct solution of the coupled differential 
equations of motion of the blade is such atr^ap"* 
proach, permitting optimization studies at low 
cost. This paper presents a model of a nonunl-* 
form, tapered, twisted cantilever wind turbine 
blade and a method of solutlon- 

For the purpose of demonstrating the method, only 
the in-plane and out-of-plane bending modes are 
considered since the torsional modes occur at 
fraquencies much higher than the bending modes 
because of the high torsional rigidity of the 
blade* The coupled differential equations of 
motion were transformed into a set of first order 
equations and solved with Runge— Kutta numerical 
Integrators. The turbine blade under study has 
tuajor stiffness di cont inul t ies . The blade is 
therefore considered as if made of adjoining 
segments, each one having a varying stiffness. 

The first and second derivatives of the stiffness 
curves evidently have to be considered. The con- 
tinuity of the shear forces and moments was im- 
posed between each segments of the blade. With 
these variations of the Runge-Kutta method it is 
possible to obtain the resonant frequencies and 
the normalized distributions of displacement, 
bending moment and shear force for the first six 
beading modes. The coriolis forces and the 
tension force due to centrifugal loading are in- 
cluded. The modes of vibration are computed for 
a regime of rotational speed. The normal and 
tangential aerodynamic loading at each section of 
the blade could be included as extra terms in the 
differential equations; however, this paper con- 
siders only a rotor turning without the aerody- 


namic forces. These will be included in further 
studies dealing with the optimization of small 
capacity wind turbines. 

An important extension to this model consists In 
the formulation of the polar moment of inertia of 
the rotor in terms of a couple at a frequency U3 
applied to the hub and the resultant angular ac- 
celeration of the rotor. A direct method using 
the shear forces _and moments produced by each 
blade at the hub and an integral from using the 
in-plane displacements (mode shape) along the 
blades are formulated and the numerical results 
compared. The polar moment of inertia is then 
obtained as a function of frequency and rota- 
tional speed. 


THE EQUATIONS OF MOTION AND THE TRANSFORMATION 
METHOD 

Lets consider a turbine blade turning outside the 
aerodynamic and gravity field with its main axis 
perpendicular to the rotation axis mounted as a 
cantilever into a rigid hub. The tension in the 
blade due to centrifugal loading is included. If 
Che blade elongation is assumed to be small re- 
lative to the transversal displacements. It can 
be demonstrated that the coriolis forces become 
negligeable. When only the In— plane and out-of- 
plane bending motions are considered, the coupled 
differential equations become: 



These equations are a subset of the equations of 
Houbolt and Brooks (ref. 1) and an extension of 
the equations of Canergie and Dawson (ref. 2) for 
a twisted blade. In these equations, the tension 
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T at any section of the blade is independent of 
the vibration frequency but proportional to 


(5d) 

and represented by: 

T 


(5e) 

2 ^ 

T = /* m(z)z 3 z ( 2 ) 

*z 

n-h 

(5f) 

Taking the first and second derivatives of the 
terms in braket while considering the moments of 

Y^ .Y 3 

(5g) 

inertia Ixx> ^yy ixy variable along Che z 

Yi - Y, 

(5h) 

axis, the system of equation ( 1 ) becomes: 
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They can be conveniently represented by: 
Y- - (z, Yj. Yg) 

i* 1, 2, •••> 8 
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where ( )’ and ( )"... Indicate the first, sec- 
ond, derivatives with respect to z. It is 
seen from these equations that the in-plane and 
out-of-plane bending are coupled through 1 ^ and 
its first derivative I^y 


The method of transformation of variables, first 
proposed by Canergie and Dawson (ref- 2) is gen- 
eralized by applying it at any frequency of vi- 
bration and not exclusively to find the resonant 
frequency- The two fourth order equations are 
transformed into eight first order equations by 
the following substitutions: 



Y 3 - v" - V ( 4 ) 

Y4 « U" Yg = U 


Applying the appropriate boundary conditions, it 
is then possible to solve equation ( 6 ) with the 
use of Runge-Kutta numerical integrators. 


The displacements and their first, second and 
third derivatives are then obtained- From these, 
one can compute the distributions of the shear 
forces and moments by the usual relations: 
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THE BOUNDARY CONDITIONS 


By substitution and differentiation the following 
eight first order equations are obtained: 


El Y’ + El Y^ 
xy 1 yy 2 


+21’ Y* + X" Y. 
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Consider the turbine blade shown in figure 1. 
The longitudinal axis of the blade is the z axis 
and the rotor turns in the y-z plane. The x axis 
is the axis of rotation. The blade is mounted in 
a rigid hub at 6 % of its span ( jO and has a total 
length of 4*95 meters (L). 


Because the blade is rigidly mounted Che boundary 
conditions are: 


EI^ y: + El Y. 
XX 1 xy 2 


2 = [2 

^ ’;;y "4 Yj . Y3J 


xy 2 


+ T Y 3 + T’Y^ + m 


y* * y 

Y3 Yj 


(5b) 

(5c) 


u = V”U*=v’=“0 at z ^ Z 

tha^ is Yi(X) - 0 for 1 - 5. 6 , 7, 8 ( 8 ) 

At the free end the conditions are; 

u" = v" - u"’ v"’ * 0 at z = L 

since the moments and shear forces are zero; 
that is Yi (L) « 0 for 1 » I, 2, 3, 4 (9) 
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Figure I - Geometry of IREQ turbine blade 


These boundary conditions apply to a stationary 
blade, l^en the hub is allowed to rotate around 
the X axis the conditions at the free end remain 
unchanged but the displacement v and its slope v’ 
are different from zero at z = However, u and 
u* remain zero since the hub is assumed rigid. 
Ther^ore, if the rotor Is allowed to spin and 
oscillate in the plane of rotation the boundary 
conditions at z = ^ become: 


= u = u’ =0 

- y * v(^) (10) 

~ - 1. ^ 

r' ThisTTast condi ti on Is the consequence of the 
“ rigid hub being displaced by an amount v(il) at 


THE METHOD OF SOLUTION 

The loethod of solving equation ( 5 ) is as follows: 

1) A value is selected for the frequency 
o>* 

2 ) The four known boundary conditions 
(equation 7 ) at the root are set and 
the four unknown conditions are given 
arbitrary values namely: 


Yj = C„. = Y3 = = 0 ( 11 ) 

3 ) From these eight boundary conditions at 
the root a solution is obtained with 
the use of fourth order Runge-Kutta in- 
Cegrators. Eight values, to Y3, are 
obtained for the free end: 

= h,i 

(r = 1) (12) 


4 ) Step ( 3 ) is repeated by successively 
setting each of the unknown boundary 
condition to the arbi t fafy“value Co- 
in this way four sets of starting 
boundary conditions are obtained for 
z = 1 : 
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5 ) 

The ! 

solution 

of equation (6) is a 


combination of these four solutions. 
Name ly , 


r=*l * 

(1-1,2, 8) 


( 15 ) 


However, the known boundary conditions at the 
free end are for a cantilever blade: 


Y. (L) =0 (i = 1 , 2 , 3 , 4 ) 


( 16 ) 


where the subscript r - I indicates the since the shear forces and moments must be zero, 
solution with the first group of bound- The right hand side of equation ( 15 ) can then be 
ary conditions; partitioned: 


179 


4 

X; F. - 0 (i = 1,2, 3, 4) (17) 

Si ^ 

6 ) A non-trivial solution is possible if 
the determinant of the coefficients 

equation (17) is equal tO zero 

>0 (i » 1 , 2 , 3, 4) 

(r = 1 , 2 , 3, 4) (18) 

Therefore, the above steps are repeated with in- 
creased values of the vibration frequency (x) until 
equation (18) is satisfied- That 0 ) then corre- 
sponds to a resonant frequency. 

7) Having found the resonant frequency 
a^ is set to 1 and 32 , 33 and 34 are 
computed giving the four unknown 
boundary conditions at the anchor 
point - 

8 ) The solution is repeated once more with 
the following initial conditions; 



Y 3 > 83 =. 0 (19) 

In addition to the displacement u & v the shear 
forces and moments are computed at each blade 
station with equations (7a) to (7d). 

The method described above was first used by 
Canergie and Dawson (ref. 2) to find the natural 
frequencies of a straight constant section blade. 
Its application here Is extended to twisted, 
tapered blade having discontinuity of rigidity. 



TURBINE BLADE CHARACTERISTICS 

The turbine blade used on the 40 kW, 10 meters 
IREQ wind turbine is a twisted, tapered composite 
blade made principally of steel and fiberglass. 
Figure i shows its construction. Its assymetric 
aerodynamic profile is NACA 4415. The chord is 
44-45 cm at the root and 10.92 cm at the free end 
with a thickness varying from 7.11 cm to 1.78 cm. 
The twist angle 3 goes from 47.2 degrees to 2 
degrees at the tip- 

The principal moments of inertia Ixxg ^YYs 
for a group of typical blade sections were com- 
puted from an engineering drawing of the blade 
and were transformed into the blade principal 
axis (in-plane and out-of-plane) by the usual 
relations: 


I 

XX 


I sin^a 

yy. 


I » I cos^ 6 I sin^B 

yy yyg 


8LRDE CHflRflCTERISTICS 



Figure 2 - IREQ turbine blade characteristics 



The Ixxg and lyyg for each element of the blade 
section are equivalent moment of inertia based on 
the same reference modulus of elasticity E. The 
blade actual geometry was found to be signifi- 
cantly different From the drawing geometry. The 
moment of inertia Ixxg lyy^ are then cor- 
rected to account for these manufacturing inac- 
curacies. The values used as input to the modal 
analysis program are shown in figure 2. The root 
sections of the blade is approximated by linear 
distributions while the all fiberglass section 
from the end of the steel Insert to the blade tip 
was approximated by a polynomial computed from 
twelve input data points. The first and second 
derivatives, of Ixx* ^yy I were numerically 
computed from the above distributions. 

In order to take care of the discontinuities of 
stiffness, the blade is divided into three sec- 
tions, the first from the anchor point to the 
blade root section, the second up to the end of 
the steel insert and the third to the tip of the 
blade. 

The numerical integration is done from the anchor 
point to the tip of the blade in a continuous 
manner except that the values of four of the 
eight state variables (Yi to Y 4 ) are varied in a 
stepwise manner at the two major discontinuities 
of the blade. This is done because the state 
variables being integrated at each segment of the 
blade are the numerical derivatives of the dis- 
placements u and V, namely u’ , v* , u", v", u"* and 
v"* and not the forces and moments in the blade. 
The physical quantities that must be continuous 
are the displacements, the slopes, the moments 
and the shear forces and the tension, not the 
derivatives. Therefore the continuity of the 
boundary conditions become; 
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The Indices - and + represent the sections imme- 
diately to the left and immediately to the right 
of the discontinui ty • Using equation (7a) to 
(7d) and equation (20)> it is possible to find 
new values (the + values) for the derivatives u", 
v**j u'% and v'” that will ensure the continuities 
of the forces and moments, across the dis- 
continuity. Performing the appropriate algebra 
on equations (7) and (20), we get: 




it is clear that the shear forces and bending 
moments computed with this method are not exact 
in the Imraediate region of the discontinuities. 
However the distribution should not be affected 
in regions farther from the discontinuities. 

FORMULATION OF THE POLAR MOMENT OF INERTIA 

The polar moment of inertia J of an horizontal 
axis wind turbine rotor is required for the 
analysis of the dynamic torsional stability of 
the drive train. J varies with the vibration 
frequency U) and the rotational speed * The 



Figure 3 - Internal forces and moments at the 
blade anchor point. 


variations of J are mainly caused by the trans- 
verse vibration of the blades. The polar moment 
of inertia of the rotor, excluding the hub, is 
defined at a frequency o) by; 



where Cv Is the amplitude of the coupl^ applied 
by the hub at the frequency w and -0)^0 is the 
angular acceleration of the hub. For a sym- 
metrical three bladed rotor, J will be three 
times the polar moment of inertia of one blade 
computed with respect to the rotor hub. 

When a blade vibrates at a frequency u) , an in- 
ternal shear force and moment appear at the 
anchor point as shown in figure 3. The shear 
force and moraent produced by the out— of— plane 
vibrations are reacted upon by the rigid hub and 
do not appear in the formulation of J. However, 
the hub is free to rotate around its axis and the 
in-plane vibrations will be reacted by the rotor 
hub In the form of a couple. The sign conven- 
tion shown in figure 4 is introduced such that 
F^ - - iM /^Z. The external couple , applied 
by the hub is then: 

C - - F (^) ft - M (ft) (23) 

V v v 

The tension T, being purely radial at the hub, 
does not produce any couple. The external couple 
expressed by equation (23) is then Introduced in 
equation (22) to give: 



Figure 4 - Positive directions of shear forces 
and moments 
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Since T is radial at F (£) and M (it) 

= [- k ^ ^ 

M (i) « El v" + El u'* 

V XX xy 

L ^ J z = i 


become : 
(25) 

= Z 


(26) 


solved in the manner described above except that 
a value is chosen for to and the arbitrary value 
C 4 used as initial condition for (v"') is 
varied until the determinant (equation 18) 
becomes zero. When a solution is obtained for 
that to, equations (30) and (31) are used to com- 
pute J. This procedure must be repeated for each 
value of to with the following boundary conditions 
at z Z I 


If the blade section is uniform and symmetric at 
the anchor point (a circular section for the IREQ 
blade) we have; 

r (Z) = a) = I (^) = 0 (27) 

XX xy xy ' ' 

and (25) and (26) become: 
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These boundary conditions were explained earlier. 


giving 


J 


jV(i) - iw-a) 
0 )^ 


(30) 


In this last equation, the polar moment of iner- 
tia of one blade Is expressed in terms of the 
forces and moments at the rotor hub. 

Another representation of the polar moment of 
inertia is possible if one considers the in-plane 
displacements of the blade as it vibrates at a 
frequency u). By using integration by parts on 
equations (1) and (7) and applying the 
appropriate boundary conditions, It can be 
demonstrated that 


T ^ r ra V z dz 

" v(^) J 

Z 


(31) 


At a very low frequency, the blade is not de- 
formed and moves as a rigid body; v(z) becomes a 
straight line 

v(z) = v{Z) z/Z (32) 

and (31) takes the well known form of the static 
moment of inertia: 


Jo 



dz 


(33) 


Both equations (30) and (31) can be used to 
compute the polar moment of inertia but the 
integral formulation is inherently more exact 
from the numerical point of view because it only 
uses the blade in-plane displacements while 
equation (30) uses in addition the second and 
third derivatives of these displacements at the 
hub anchor point. 


NUMERICAL RESULTS FOR THE MODES OF VIBRATION 

The natural modes of vibration for the IREQ HAWT 
blade have been computed for the following 
conditions; 

1) The first six modes of a stationary 
cantilever blade. 

2) The first six modes of a cantilever 
blade at 100, 200 and 300 RPH. 


FIRST MODE 

rneouENcr. 2. B 72 za Hcnri 





Figure 5 - First out-of-plane mode 


In order to compute the polar moment of inertia, 
the equations of motion of the blade must be 
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SECOND MODE 

FNCauCNCT- a. 3672S ttCflTZ 



'CO O.i 0.2 0.3 0.*i 0.3 0.9 Q.7 0.8 0.9 l.O 



'C.O O.L 0.2 o '. 3 0.9 0.3 0.6 0.? 0.8 0.9 l.O 



Figure 6 - Second out-of-plane mode 


THIRD MODE 



’0 0 0,1 0.2 0.3 0.9 0.5 O.S 0.7 0.8 0,9 1.0 



'0,0 0.1 0.2 0.3 0.9 O.S 0.6 0.7 0.8 0.9 l.O 
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Figure 7 


First In-piane mode 


FICE END Of TuneiNE OLflOC ftt£ END Of lUttlNE 


FOURTH MODE 

FReOUCNCT- 21.H1S38 HCNTZ 



’0.0 0.1 0^2 O’S 0.9 0.5 0.6 0.7 0.8 0.9 l.O 



Figure 8 - Third out-of-plane mode 


For each of these cases, the following has been 
obtained : 

1) The resonant frequencies (poles). The 
zeros are also available from the polar 
moment of inertia curves. 

2) The normalized In-plane and out-of- 
plane displacement curves or mode 
shape. 

3) The normalized in-plane and out-of- 
plane shear force and bending moment 
distribution curves. 

The results presented here are valid only for a 
blade mounted perpendicularly in a rigid hub. 
Only the coupled in-plane and out-of-plane 
bending modes are considered, torsion being 
neglected. One blade was tested experimentally 
in the laboratory for the stationary case only. 
The analytical and experimental frequencies are 
compared in table I . 


^^0DES 

ANALYTICAL 

EXPERIMENTAL 

1 

2.872 

2.80 

2 

8.387 

8.00 

3 

iO.627 

10.99 

4 

21.415 

18.66 

5 

31.384 

27.39 

6 

37.474 

30.77 


Table I. Natural resonant frequencies in Hertz 
for a stationary blade. 




It can be seen that the natural frequencies are 
in close agreement for the first few modes- The 
larger discrepancies for the higher modes are 
believed to be caused by some uncertainty in the 
construction causing local variations of mass and 
stiffness which would affect mostly the higher 
modes- Also, the fact that the blade support was 
not perfectly rigid, means that the observed fre- 
quencies would be lower than the one predicted 
under the cantilever assumption- Finally, the 
excitation of the blade was done with an electro- 
magnetic exciter which requires some attachment 
hardware at the tip of the blade. This addition 
of mass at the tip, also tends to lower the f re- 
queue les- 

The normalized mode shapes, shear forces and 
bending moments distribution curves are shown in 
figures 5 to 10 for a stationary blade* Coupling 
between the in-plane and out-of-plane modes is 
evident from the figures. It should be noted 
that the blade stiffness discontinuities at blade 
station O.IO and 0.22 do not affect the conti- 
nuity of the distribution of the bending moment 
and of the shear force. 

The same computations were done for a rotating 


blade 

at 

100 , 

zoo and 300 

RPM. 

Table 

II 

shows 

MODE 

RPM =* 0 

RPM 

* 100 

RPM 

= 200 

RPM 

= 300 

1 

2 , 

872 

3. 

728 

5. 

.421 

7, 

,183 

2 

8 . 

387 

8 . 

649 

9, 

.099 

9. 

,594 

3 

10 . 

626 

11 . 

228 

13. 

.031 

15, 

.703 

4 

21 , 

.415 

22 . 

143 

24. 

.161 

27, 

,117 

5 

31 

.384 

31. 

593 

32. 

.206 

33, 

,143 

6 

37. 

474 

38. 

206 

40. 

.246 

43, 

.486 


Table II. Resonant frequencies of a rotating 
blade. 

the resonance frequencies obtained. The effect 
of the rotational speed on the resonance fre- 
quencies for the three first modes is shown in 
figure 11. The agreement with some experimental 
results is good. 


NUMERICAL RESULTS FOR THE POLAR MOMENT OF INERTIA 

The polar moment of inertia of a three bladed 
rotor has been computed for vibration frequencies 
up to 70 radians/sec and rotational speed up to 
200 RPM, The values of Jp -j. obtained with the 
formula using the shear force and moment (equa- 
tion 30) and the values Jf obtained with the 
Integral formula (equation 31) give comparable 
results. However, as mentioned earlier the 
values of Jf are implicitly more accurate than 
the values of Jp 4 . This fact is demonstrated 
numerically by observing that at very low fre- 
quencies, the value of Jj remains constant for 
all RPM used while the value of Jp + shows 
small variations for each rotational speed 
considered. 


FIFTH MODE 


^REOUCHCT. 31 . 38390 Mtnri 



'a.O 0^1 O.Z 0.3 G.(« 0.5 0.8 0.7 0 8 O.f t.O 



'0.0 O.t 0.2 0.3 0.>4 0'5 0.8 0.7 0-8 0.9 l.Q 



Figure 9 - Second in-plane mode 


SIXTH MODE 

FftCQUEMCr. 37.47387 ME8TX 



'0.0 0.1 0.2 0.3 0.4 0.5 0 6 0.7 3.9 0.9 1.0 



Figure 10 - Fourth out-of-plane mode 
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ANALyilCAL AND EXPEfilfOTAL 


0 Ahalytical 
• Experimcntal 



Figure ll - Effect of rotational speed on vi- 
bration frequencies 


POLAR MOMENT OF INERTIA OF ROTOR 



Figure 12 - Variation of polar moment of inertia 
In fonction of vibration frequency 
and rotational speed 


The variation of the polar moment of inertia with 
the vibration frequency and rotational speed is 
presented in figure 12. Only the values of J; 
are shown. The poles (resonant frequencies) and 
the zeros can be seen on this graph. The nu- 
merical value of the zeros are presented in 
table III. 


ZEROS 

RPM = 0 

RPM = 100 

RPM “ 200 

1 

2.916 

3.789 

5.588 

2 

9.912 

10.689 

— 


Table III. Zeros of the rotor in Hertz 

It can be seen that the value of J is relatively 
constant at low frequency and comparable (less 
than 0.47Z at to * 0.5 rad/sec) to the static 
value Jq. (The static value is 807.29 newton- 
laeter-sec ) 

If damping had been included in the equations the 
extreme variations of J at a pole would be re- 
duced, especially when a zero is very close to a 
pole, as is the case for the first mode. With 
damping the p6le-zero doub^let would produce only 
a small variation in J, its importance depending 
on the separation between the pole and the zero. 

CONCLUDING REMARKS 

A mathematical model and its method of solution 
have been presented for a tapered twisted, can- 
tilever wind turbine blade with discontinuous 
stiffness. The two fourth order differential 
equations representing the in-plane and the out- 
of-plane motion of the blade have been trans- 
formed Into eight first order equations and 
solved with Runge-Kutta integrators. The blade 
discontinuities have been approximated by im- 
posing the continuity of displacements, slopes, 
bending moments and shear forces. The centrifu- 
gal force is included in the model; the coriolis 
force was found to be negligeable. The polar 
moment of inertia of a three bladed rotor is for- 
mulated considering either the in-plane bending 
moment and shear force at the anchor point or _ the 
integral of the in-plane displacements for vi- 
brating, rotating blades. 

It has been demonstrated that the method is 
sufficient to compute the natural frequencies and 
mode shapes of a stationary or rotating wind tur- 
bine blade with large discontinuities in stiff- 
ness. The normalized distributions of bending 
moment and shear force are also computed. The 
polar moment of inertia has been computed as a 
function of frequency and rotational speed. Good 
agreement with experimental frequencies has been 
observed. 

The computer program can be used efficiently for 
the structural optimisation of the blades of 
horizontal axis wind turbine. The computer time 
and memory requirements are relatively small (ap- 
proximately 20 sec and 200 K with an IBM 370, for 
each mode) so that parametric studies are pos- 
sible. 
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NOTE: 


NOMENCLATURE 


Numerical results of vibration frequencies and 
raode shapes for discontinuous turbine blades 
published by Lang and Nemat - Nasser (ref. 3) 
became known to us just recently, after the 
analysis presented here was completed. The 
accuracy of the method proposed here will be 
compared later with the results of reference 3. 
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Coefficients in numerical solu- 
tion 

Arbitrary constants used in nu- 
merical solution 

Couple applied by hub at anchor 
point of the blade 
Young's modulus 

Shear forces in x and y direc- 
tions 

Moment of inertia of blade about 
X axis 

Moment of inertia of blade about 
y axis 

Product of inertia 

Total length of blade 

Distance between axis of rotation 

and anchor point of the blade 

Bending moments about x and y 

axis 

Mass of blade per unit length 
Rotational speed in rev/min 
Tension force In blade 
Displacement along x axis 
Displacement along y axis 
Cartesian coordinates 
Variables of transformation 
Twist angle of blade 
Angular displacement of rotor 
Rotational speed of rotor 
Vibration frequency 
First derivative with respect to 
z 

Second derivative with respect to 
z 

Third derivative with respect to 
z 

Fourth derivative with respect to 
z 
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QUESTIONS AND ANSWERS 


G. Beaulieu 


From; W,C^ Walton 

Q: Would you agree that the root support stiffness should affect the lower modes first 

so that this is probably not the explanation for higher mode errors? 

A; True, The support stiffness could explain the lower frequency in the first few 
modesj while unknown mass and stiffness distribution and the tip mass addition 
of the electromagnetic exciter could explain the deviations for the higher modes. 

From; W.N, Sullivan 

Q; How were the experimental resonant frequencies shown measured on the turning 
rotor? 

A: directly measured from strain gages recording on strip charts. We would have pre- 

f erre^ magnetic tape recording and spectral analysis. 

From; A. Wright 

Qi Why do the boundary conditions for edgewise displacements change if hub is free to 
rotate? 

A: When the huh is free to rotate^ an in-plane displacement v(l) is present at the 

anchor point and similarly the slope of this displacement is v(l)/t. If the huh 
would be fixedy v and v* would equal zero. 

From: Y,Y. Yu 

Q; Could you elaborate further on the blade construction? 


A: Referring to Figure ly we can see the stall shaft and plate insert near the root. 
The steel is hounded to the fiberglass box which is present for the complete span 
_of the blade. The fiberglass box is hounded to the outside skin having a RACA,4415 
ptofi le , 

From: A.D. Garrad 

Q; Do you have an estimate for the damping in the blade? 

A: NOy but some stationary blade vibration tests are being done now and exponential 

decay will he measured to obtain an estimate of structural damping. Aerodynamic 
damping will not be measured. 
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